On the inner structure of confined positronium 
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We study the inner structure of a positronium confined to a void by treating it as a two-particle 
quantum system where the electron is captured by an infinite potential well and the positron feels 
only the Coulomb attraction of the electron. We calculate the ground state and the related quantities 
that govern the annihilation rate, and discuss implications of the results. 
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I. INTRODUCTION 



Shortly after the discovery of the positron (e + )i, it was 
recognized that e + may form an 'atom' (e + e _ ) which, 
from the point of view of spectroscopy, may be looked 
upon as an ultralight isotope of hydrogen^. With the 
same right, (e + c~) may be considered as an ultralight 
isotope of antihydrogen. Together with the photon, the 
Majorana neutrino, and the electrically neutral 7r-meson 



(c+c ) thus belongs to those 'particles' that are 



their own 'antiparticlcs'. In 1945 (e + e~), at that time 
still hypothetical, was named positronium 3 ., later to be 
equipped with the 'chemical' symbol Ps 4 . At about the 
same time its main properties were established^. Among 
these are that Ps is rather short-lived due to the particle- 
antiparticle annihilation and that its ground state splits 
into two substates, a singlet state ^So and a triplet state 
1 3 Sl 

In the ^So state, Ps is called parapositronium (p-Ps), 
and in the l 3 Si state, orthopositronium (o-Ps). Under 
most circumstances, the energy excess of o-Ps over p- 
Ps, A£„-p = 0.8412 x 1CT 4 eV, is negligibly small, 
hence the ratio of the production rates of o-Ps and p- 
Ps is usually equal to that of the statistical weights of 
the substates, namely 3:1. The most important dif- 
ferences between o-Ps and p-Ps result from the differ- 
ent annihilation modes of these two substates, which, in 
turn, are consequences of the conservation of energy and 
angular momentum. The self-annihilation of p-Ps re- 
sults in the generation of two 7-quanta of equal energies 
Ey = m c c 2 = 0.511 MeV, where m c is the electron mass 
and c the speed of light in vacuum, whereas the self- 
annihilation of o-Ps requires the generation of, at least, 
three 7-quanta. Owing to the smallness of the coupling 
constant of quantum electrodynamics (a rj 1/137) and 
to various statistical factors, the ratio between the 27 
annihilation rate, T2 7 = 7.99 x 10 9 s^ 1 , and the 37 an- 
nihilation rate, T3 7 = 7.04 x 10 6 s _1 , exceeds 10 3 and is 
therefore highly significant. 

The preceding statements refer to Ps in vacuum or 
in extremely dilute gases. In dense gases and in con- 
densed matter the fate of o-Ps is radically changed, since 



now the e + -e _ annihilation may (and in general does) 
involve electrons other than the Is electrons of the Ps 
atoms. This leads to 27 annihilation by the 'pick-off' 
mechanism. This mechanism may be visualized as fol- 
lows. An electron in the neighbourhood with spin direc- 
tion opposite to that of the e + 'picks off' the positron 
from a Ps atom and annihilates it by a 27 reaction. The 
pick-off reaction rate r po is determined by the probabil- 
ity density of the electrons with opposite spin at the e + 
location. This density is clearly less than that of elec- 
trons in the Is orbit of Ps, in condensed matter typically 
by about a factor of ten. Hence r po is smaller than the 
self-annihilation rate T2 7 but still large compared to T3 7 . 
This has several important consequences: 

(i) Positron annihilation in condensed matter takes 
place virtually exclusively by the generation of two 
7 quanta. 

(ii) The annihilation rate of o-Ps depends significantly 
on the environment in which the positronium atom 
is located. 

(iii) Because of (ii), the o-Ps annihilation rate may be 
used as indicator of the preferred sites of Ps in con- 
densed matter. 

(iv) In Ps-forming solids containing imperfections that 
attract Ps, the site-sensitivity of r po may be used 
to detect and monitor such imperfections. 

The pick-off annihilation rates of o-Ps fall into a range 
that allows accurate measurements of their inverses, the 
o-Ps lifetimes^. Often it is possible to perform e + life- 
time spectrometry, i.e., to distinguish different positron 
annihilation sites by the different e + lifetimes associated 
with them. In order to gain a theoretical understanding 
of the site dependence of the e + lifetimes, one has to cal- 
culate the wavefunctions of Ps atoms that arc confined 
to the various imperfections that are capable of trapping 
Ps. This has turned out to be a very formidable task. 

The electrostatic forces experienced by electrons and 
positrons are equal and opposite. In most situations the 
variation of the electrostatic potentials over the diameter 
of a Ps atom (about two Bohr radii, see below) is negligi- 
bly small compared with the positronium binding energy 
-Efree-ps ( m vacuum about 6.8 eV). Unless a Ps atoms 
enters a kind of hydrogen bond, the electrostatic forces 
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on its constituents cancel to an excellent approximation. 
Then the interaction of Ps with its condensed-matter en- 
vironment is entirely due to the exchange interaction be- 
tween the Ps-electrons and the host electrons since for 
c + implanted into matter the Pauli principle is irrele- 
vant. The fact that the ensuing repulsion of Ps by the 
host atoms is short-ranged suggests that the confinement 
of Ps to interstices or locations where atoms are missing 
may be represented by a rectangular potential well at 
these sites. 

Since even this model is still a nontrivial quantum me- 
chanical problem, Tao and Eldrup initiated a further sim- 
plified model in which the positronium is considered as a 
point particle in a spherical infinite potential well, being 
in the ground stat ofi?i? , extending the ideas of Ferrelli^ 
and Brandt et a/ii. In this model, pick-off annihila- 
tion happens in a layer ro < r < Ro near the surface 
of the spherical volume of radius Ro, the layer repre- 
senting the medium outside the hole. The width of this 
layer is a parameter fitted experimentally. Later, this 
model was refined by considering a potential wall of fi- 
nite height*^ and of unsharp shaped. In these modifi- 
cations, the pick-off annihilation of the positron occurs 
in the region Rq < r < oo, which seems to be a better 
approximation from the physical point of view and pro- 
vides a better fit to the experimental data. Nevertheless, 
as stressed by Mukherjee et ali&, this approach is still 
an oversimplification since it treats the positronium as a 
single point particle. 

In this paper, we provide an approach that avoids this 
simplification. We determine the quantum mechanical 
ground state of the positron-electron two-particle system 
with an infinite potential well that acts only on the elec- 
tron, in addition to the Coulomb attraction between the 
positron and the electron. Hence, we give account of the 
arising nontrivial internal structure of a confined positro- 
nium. 

Due to the complicated nature of the problem at hand, 
we calculate the ground state numerically, in terms of 
an expansion with respect to an appropriate set of ba- 
sis functions (truncated to contain the first 216 basis 
functions) and via a variational method, for various pore 
sizes. The probability of finding the positron outside the 
void and the overlap of the positron and the electron in- 
side the hole are also determined. We analyse and inter- 
pret the results, and discuss how further refinements of 
the model and the explanation of the observed tempera- 
ture dependence of the pick-off rate can be implemented. 

It is to be noted that a related investigation was done 
by Sommerfeld and Welker— , where the ground-state en- 
ergy and wave function of a hydrogen atom in a hole is 
calculated. In that study, the proton is considered to be 
infinitely heavy and to reside in the center of the spheri- 
cal hole, and the electron is confined to the cavity by an 
infinite potential well. We make a comparison between 
those results and our ones, and provide a simple common 
explanation of the behavior of both systems. 

II. THE MODEL 

Using the nonrelativistic framework and observing 
that, under the present conditions, spatial and spin de- 



grees of freedom are decoupled, we formulate our model 
expounded above with a normalized electron-positron 
two-particle scalar wave function >l/(r e ,rp), the Hamil- 
tonian 

n = -^(v 2 T + ) - / , 1 , (i) 

2rnc rc TpJ 47T£ |r e — r p | 

and the boundary condition 

*(rc,r p )| rc = flo =0, Vr p . (2) 

Here, r e denotes the electron co-ordinates and r p the pos- 
itron co-ordinates, Rq is the radius of the potential well to 
which the electron is confined, e is the elementary charge, 
Sq = 8.854187 • I0~ 12 As/Vm is the dielectric constant 
of the vacuum, and 2Trh~ = h is Planck's constant. 

We wish to determine the ground-state wave function 
— from now on, ^(ro, r p ) will refer to the ground state — 
together with the ground-state energy, and to calculate 
the following two ground-state related quantities. First, 
the probability 

Pont= J d 3 r p y*d 3 r c |*(r c ,r p )| 2 (3) 

r P >Ro r B <R 

that the positron will be found outside the hole. This 
quantity is a global measure of the pick-off annihilation 
rate of o-Ps. Second, the electron-positron contact pa- 
rameter 

/d 3 r c |vE<(r c ,r p = r c )| 2 

k = , (4) 

J d r e |*f r0c _ps(r c ,r p = r c )| 2 

where ^frcc-PS is the ground-state wave function of the 
free positronium, having the form 

7^f exp ( _l£ ^) (5) 

after separating the centre-of-mass motion, with a# = 
Aneoh 2 / (m G e 2 ) = 0.529 x I0~ 10 m denoting the Bohr 
radius. This quantity measures the overlap of the elecron 
and the positron, and is equal to the ratio of the self- 
annihilation rate of confined Ps to that of free Ps. 

Because of the spherical symmetry of the system, the 
ground-state wave function depends on three indepen- 
dent scalar variables only. For these we choose the dis- 
tance of the electron from the hole centre, r c = r G |, 
the electron-positron separation R = |r p — r e |, and 
the angle x between r e and r p — r , satisfying cosx = 
r ■ (r p - r c )/(r c R). 

A lengthy yet straightforward calculation shows that, 
in terms of these variables, the Hamiltonian reads 



3 



n 



+2 



2m, 
sinx 



dr 2 

d 2 



if 



2rl + R 2 



dR 2 

2_ d_ 

dRdx r e dr e 



2r c Rcosx d 2 
d\ 2 

\- R 2 cos x 



r 2 R 2 
4 d 2r 2 cos % 
RdR + 



2 cos x 



i) 2 



\R 2 sin% 



dr c dR 
2r c R d 
dx 



R 



d 2 



dr c dx 
1 



47T£n R 



r 



(6) 



The Hamiltonian @ is not separable and does not con- 
tain a small parameter on which a perturbation expan- 
sion might be based, therefore, we perform the subse- 
quent calculations numerically. 



III. THE METHOD OF COMPUTATION 

We use a variational method to compute the approxi- 
mate ground-state wave function, minimizing the energy 
(^f,7i^f) in the function space 



N 



]T c n # r 



(7) 



where the functions # n (n = 1, . . . , N) form a set of ap- 
propriately chosen orthonormal base functions. The min- 
imization is carried out by finding the eigenvector of the 
matrix of the elements Tt m n = (^imWn) with the low- 
est eigenvalue. The components of this eigenvector are 
then identified with the coefficients C n of the approxi- 
mate ground state wave function (J7J . The base functions 
are chosen in such a way that 

(i) the r -dependence of the wave function is similar to 
the ground state of a single electron in an infinitely deep 
potential well, 



t/je{r e ) = const. 



sin(7rr e /i?o) 



(8) 



(see, e.g., the work of Galindo and PascuaU^), 

(ii) the /independence resembles that of the free positro- 
nium (cf . Eq. [SJ) , and 

(iii) that the matrix elements can be computed relatively 
fast. 

To this end, we introduce the functions 



V ( i>(r e )= {r 2 c -R 2 y (i=l, 
<pf\R)= iP'^e-*/* (.7 = 1, 
Vk\x) = cos^x (fc = l, 



Ni), (9) 
N 2 ), (10) 
, N 3 ), (11) 



from which, by means of the Gram-Schmidt orthonor- 
malization procedure 1 ^, we obtain the orthonormalized 
functions ^\r c ) , ip^ 2 ' ' (R) and t/'i (x) > respectively. 
The products 



* ijk (r e ,R, X ) = ^\r c ) ■ ^'(R) ■ V-rW ( 12 ) 



(2), 



',(3)/ 



serve as the base functions \l/ n , i.e., the expansion Q is 
implemented as 



JVi N 2 N 3 
i=l j=l k=l 



(13) 



The variational method is improved by considering £ 
as an additional, nonlinear adjustable parameter, with 
respect to which the energy is also minimized (for a free 



positronium, £ = 2ob = 1.058 x 10 10 m 



For the first functions ipi(r e ), ip\ z> (R) and Vi^Cx)) 
the orthonormalization procedure means only a normal- 
ization. It is plausible to expect that the product 

tfm = const. ■ #(r e ) ■ ^(R) ■ V f\ X ) (14) 

proves to be the dominant basis function in the ground 
state JTJJ. 

The matrix components Timn{Ro,C) were calculated 
analytically by the software Maple V— and subsequently 
exported to the format of programming language C, and 
the numerical solution of the eigenproblem, as a func- 
tion of i?o and £, was performed by a C program using 
routines from Numerical Recipes^. For a fixed Rq, the 
lowest eigenvalue was minimized as a function of £, and 
the eigenvector corresponding to this eigenvalue was de- 
termined. The lowest eigenvalue gave the approximate 
energy of the ground state, while the components of the 
eigenvector provided the coefficients Cijk of the approx- 
imate ground-state function (cf . Eq. I13|) . Finally, using 
the resulting C^k values, p ou t and k were calculated again 
by Maple V. The computations took approximately two 
months on a workstation. 



(2) 



IV. RESULTS 

It was found that a base-function set with N\ = N 2 = 
N3 = 2 sufficed for determining the energy, Ni = N 2 — 
N 3 = 3 for p out , and TYi = 7V 2 = yV 3 = 4 for k - 
in the sense that any further increase of the number of 
base functions produced only a negligible change in the 
values of E, p OVL t, or k. The results presented below were 
calculated using the base-function set Ni = N 2 = N 3 — 
6. The calculations were performed for 33 different values 
of the hole size in the range ae < Ro < 98ob- 

The resulting ground-state energy E is shown in Fig. ^ 
as the function of the hole radius. For large Rq values, E 
approaches the ground-state energy -Ef rcc -ps = — 6.80eV 
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FIG. 1: The calculated ground-state energy of the confined 
positronium, E^ num \ as the function of the hole radius. The 
data points are fitted by the function given in Eq. JUiJ (solid 
line) . The dotted line shows the ground-state energy of a sin- 
gle electron confined to a hole, E c = Tv 2 h 2 / '(2m c _Ro) > and the 
dashed line denotes the ground-state energy of a hydrogen 
atom in a hole. The ground-state energy of the free positron- 
ium, Sfvee-ps, and the ground-state energy of the free hydro- 
gen atom, 2£f rcc _ps, are also indicated. Inset: the difference 
between the ground-state energies of the confined and free 
positronium on a log-log scale, together with the fit 115H . 



20 40 60 80 100 

R /a B 

FIG. 2: The weights of the most dominant base functions in 
the ground state, as functions of the hole radius. (Concerning 
the coefficients themselves, Cm and C121 are positive, and 
C211, C112 and C122 are negative, for all Ro.) The weights of 
■01 ( r e) and %l)^\r c ) in the ground state of a single electron 
in a hole ||SJ — approximately 0.983 and 0.017, respectively 
— are also displayed. Inset: the parameter £, as the function 
of the hole radius. The value of £ for the free positronium is 
2cib. 



of the free positronium. The dependence of the energy 
on Rq is well described by the function 



E = IE 



free- PS I 



8.75 (*> 



ii.78 (^y 



-2.53 



1 



The error of the numerical parameters in this fitted func- 
tion is less then 3.7%. For comparison, Fig. displays 
the ground-state energy of a single electron confined to a 
hole, E e = ir 2 h /(2m e R^) (see, e.g., the work of Galindo 
and PascuaU^), and the ground-state energy of a hydro- 
gen atom in a hole with its proton in the centre (repro- 
duced from the work of Sommerfeld and Welkeui^) as 
well. 

The determined ground state \& of the positronium in 
a pore is characterized by the coefficients Cijk and £, Fig. 
121 presents the most dominant coefficients and £, as the 
function of the hole radius. The coefficients are displayed 
in their squared form, ICyfcl 2 , which give the physical 
weights of the corresponding base functions ^ijk in the 
ground state 'J? . All coefficients C^k were found real, as 
was to be expected. 

The calculated independence of the probability p ou t 
that the positron can be found outside the void is shown 
in Fig. |3J The data may be fitted by the function 



1 + 0.018 



Pout 



(a B ) 



0.34(f)" +0.004 (ny 



(16) 
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FIG. 3: The probability p out that the positron of the confined 
positronium can be found outside the hole, as the function of 
the hole radius. The data points are fitted by the function 
given in Eq. 1161 . 



Similarly, the results for the density parameter k (cf. 
Fig. may be described by the function 



0.4i«!i 



0.024 



(a B ) 



3.22 



(17) 



with an uncertainty of the coefficients of less than 4% 



with uncertaities of the fit parameters of less then 2% . 
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FIG. 4: The density parameter k of the confined positronium 
as the function of the hole radius. The data points are fitted 
by the function given in Eq. Q17[l. 



V. ANALYSIS OF THE RESULTS 

The first observation to make is that, as expected, for 
large hole radii, the system tends to the free positronium. 
This can be seen on each the quantities E, £, p ut and 
k. Also, as expected, it is among the base functions 
^ijk whose presence is the strongest in the determined 
wave function, and this dominance increases with the size 
of the hole to almost 100% . 

To investigate the ground state more closely, let us 
study its dependence on the variables r c , R and x- This 
can be done by calculating the weights of the base func- 
tion components ip^\r e ), tpj(R) and tp^ (x) m tne 
wave function. These weights are defined as 

N 2 N 3 

^ (1) = EEfei 2 ( 18 ) 

3=1 k=l 

for a i>l (r e ), and analogously for a ^\r) and a 

i^u\x)- The calculated weights are shown on Fig. 
as functions of the hole radius. 

Concerning the r c -dependence, we can observe that the 
hts W t {1) are almost independent of i?o, and that 

they are very close to the weights of i/}^(r e )a in the 
ground state of a single electron in a hole [cf. Eq. ©], 

W-^- 0.983, M^ 93555 ' 1 "-* 2 ^ 0.017 

(19) 

(the higher weights | C3 1 2 , | C4 1 2 , ... are negligible so it is 
enough to consider the first two ones). Therefore, the 
r e -dependent part of ^ is close to the wave function of 
the confined electron, even for small values of i?0i an d 
with a better and better accuracy as Rq increases. 

The i?-dependence is more influenced by Rq- In a small 
void, the components higher than the first one contribute 



FIG. 5: The weights of the base functions ip[^ (r c ), ip^ ( r e), 

ip{ 2) (R), ip { 2 2) (R), and Vl 3) (x), ^2 3) (x) in the ground state 
of a confined positronium. For each of the three variables, 
the presence of the higher order base functions is small, be- 
yond the fourth order they can be practically neglected even 
for small hole sizes. The dashed lines show the weights of 
functions ^J^(r e ) and ip^ ( r c) in the ground state of a single 
confined electron (cf. Fig. |5J . 



up to 20% and this amount decreases to zero as the hole 
size is increased. For large Rq, only the first component 
remains, and in the meantime its parameter £ approaches 
the value of the free positronium. The x-dependence con- 
tains a 10% amount of the higher base functions for small 
Rq, which decreases to zero for increasing hole radius. 
Since ipi^(x) = const. , this means that, for larger pore 
sizes, the ground state becomes ^-independent. 

Turning to the ground-state energy, a simple interpre- 
tation of its found i?o _ dcpcndcnce can be that the pres- 
ence of the finite hole increases the energy of the electron 
by the amount of E c — n 2 h 2 /(2m e R l ), similarly as it 
does with a single electron. Fig. shows that this inter- 
pretation is fairly good. In Appendix A, we show that, 
with an only somewhat more involved physical argument, 
it is possible to give an even more precise approximation 
to the energy as the function of i?o- 

The probability that the positron can be found out- 
side the hole, p ou t, is one for holes much smaller than 
ae, since the characteristic space region needed for the 
positron is of the order of ob, like in the free positron- 
ium (see Appendix A as well on this subject). On the 
other side, for large hole sizes, p ou t falls down to zero as 
i?o 3 (cf. Eq.nSJ. This latt er can also be understood, an 
explanation for this is presented in Appendix B. 

The density parameter k also behaves as expected. For 
Rq large, it is 1 since the positronium is almost free. By 
gradually squeezing the hole, the allowed space for the 
electron becomes smaller, it can no more go so far away 
from the positron as in a free positronium. This results 
in some increasing in the overlap-measuring k . We can 
guess what happens when Rq becomes much smaller than 
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FIG. 6: The ground-state energy of the confined positronium 
as a function of the hole radius, together with two approxima- 
tions to it, _E froc _ps+7r 2 ?l 2 /(2mc-Ro) and Ea. lA9l The ground- 
state energy of the confined hydrogen atom as a function of 
the hole radius, together with its approximation Eq. IA11I is 
also shown. 



the Bohr radius: The electron is practically confined at 
the origin, the wave function of the positron becomes 
similar to what it is like in the antihydrogen atom, and, 
therefore, the independent part of \& becomes 



1 



R 

exp I 

a B 



(20) 



The k value corresponding to this wave function is 8 [cf. 
Eqs. and J5J] so the increase that can be observed 
in Fig. 0] is expected to slow down for Ro < ae and to 
converge to the value 8 reached at Ro — 0. 

We mention that a similar behavior of n is expected 
if one replaces the infinite potential well by some other 
confining potential. Limiting the electron to a region is 
plausible to cause a kind of squeezing the two particles 
onto each other, irrespective of the concrete form of the 
confining mechanism. 



VI. DISCUSSION 

In the framework of our model, the lifetime of an or- 
thopositronium is determined by the pick-off annihila- 
tion rate, proportional to p ut, and the intrinsic 37 decay 
rate of the orthopositronium, which is considered to be 
its vacuum value multiplied by k. The angular correla- 
tion of the two emitted photons of the pick-off process 
can be obtained by generalizing the Fourier approach of 
Mukherjee and coworkers^ to the total momentum of 
the present two-particle system. 

The infinite potential well, chosen in our calculations 
presented here, can be a good approximation for the con- 
fining mechanism for the electron in the pore for some 
materials and less good for others. The same variational 



method can be applied for other potential shapes, modi- 
fying only the electronic base functions, in the appropri- 
ate way. Naturally, finding and using the corresponding 
electronic base functions for various potential shapes may 
make the calculations hard to carry out in practice. In 
this respect, it is worth mentioning a feasible alternative 
approach for this purpose. It is known from quantum 
mechanics that the boundary condition of requiring the 
vanishing of the wave function *F, Eq. [3 is not the only 
possibility for confining a particle to a finite region 19 . 
There exists actually a one-parameter family of quan- 
tum mechanically allowed boundary conditions, of the 
form ^ + L d^> / dr c = , that is, the vanishing of a com- 
bination of the wave function and its derivative, where L 
is an arbitrary (real) length scale parameter. The well- 
known and most frequently used case, the vanishing of 
the wave function itself is simply the special case L = 0. 
Instead of choosing different potential shapes it is sim- 
pler to choose only a different value of L in the boundary 
condition. The corresponding electronic base functions 
can still be chosen as polynomials, fulfilling now the new 
boundary condition. One can expect that the precise 
shape of the confining potential is not so important in 
practice. Bearing in mind that these boundary condi- 
tions express nothing but that the quantum probability 
cannot flow outside from the holei^, and that any con- 
fining potential shape wishes to ensure the same — al- 
though not completely strictly — , we can guess that all 
confining potentials will, in their effect, resemble one or 
another boundary condition (potential well) with an ap- 
propriate L. Therefore, the possibilities provided by the 
general potential shapes may be fairly well represented 
with the one-parameter family of sharp potential wells. 

In this paper, we have calculated the ground state only. 
To give account of the observed temperature dependence 
of the pick-off decay rate, it would be useful to deter- 
mine the excited states, too. Indeed, as demonstrated by 
Goworek and coworkers^ and Gidley and coworkers^! 
for some variants of the Tao-Eldrup model, this temper- 
ature dependence can be attributed to the fact that it is 
also possible for the positronium to be in one of its ex- 
cited states, approximately with the thermal equilibrial 
probability e~ E / kT . Explaining the effect in our frame- 
work, the excited states possess different p ou t and, con- 
sequently, a different annihilation rate. Calculating the 
excited states is also possible in the variational method, 
after suitable modifications. 

In the end we mention that it would be interesting to 
perform a calculation similar to the one presented here 
for systems when the two particles are of different mass, 
e.g., when the positron is replaced with a proton or a 
fi + particle, as an interpolation between the positron 
and the infinitely heavy limit discussed by Sommerfeld 
and Welker. Unfortunately, our results cannot be simply 
'renormalized' to be applicable to such situations. The 
reason for this is that, for different masses, the Hamilto- 
nian has a more general form than Eq. © has, and these 
differences are not only some easily rescalable numerical 
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factors but mean a more general and complicated struc- 
ture of the Hamiltonian. However, the considerations of 
Section 4 and the Appendices may more easily be ad- 
justed to the case of different masses, and in this way at 
least some approximate information could be obtained 
for those systems. 
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APPENDIX A: ESTIMATING THE GROUND 
STATE ENERGY 

A closer, yet not involved interpretation of the found 
i?o-dependence of the energy of the positronium in a hole 
is motivated by the simple derivation of the ground-state 
energy of the free hydrogen atom. There, based on the 
uncertainty relation, one approximates the average mo- 
mentum of the electron by p = h/r, where r is the size 
of the average space region "run" by the electron. Thus 
the energy of the electron is 



V 



1 



1 



1 



2m 47T£o f 2m c r 2 r 



(Al) 



Minimizing this expression in r one reaches just the cor- 
rect ground-state energy, and the corresponding r is also 
nothing else than the Bohr radius a B . 

As a next step, let us consider the free positronium. 
There, the average distance between the electron and the 
positron is 2a B (cf. Eq. EJ. Neglecting the motion of the 
center-of-mass, the average region run by the effective 
particle of reduced mass m e /2 belonging to the relative 
motion is then 2a B . The energy is, therefore, considered 
as 

% 2 1 c 2 1 

(A2) 



2(m e /2) (2o B ) 2 47T£ 2a B ' 

which proves to be just the correct energy value -Ef rC c-PS- 
If we return to the 'two particles'-picture, the kinetic 
energy in Eq. i|A2(l is shared by the two particles so each 
possesses the half of it, 

1 % 2 1 £ 2 i 

(A3) 



2 2(m e /2) (2a B ) 2 2m c (2a B ) 2 ' 

Using this in the reverse way, we find that the character- 
istic space region for both particles is 2a B . 

Now, when we confine the electron in a hole, its char- 
acteristic region will be influenced by the hole as well. To 
this end, let us determine the characteristic space size of 
a single electron in a hole with the same logic as before, 
from its known energy: 



E e . = 



r 2 h 2 



l 



2m c R 2 2m c {Ro/Trf 



(A4) 



"implies" that the corresponding characteristic distance 
is Ro/nM 

When the electron is both confined to a hole and at- 
tracted by a positron, its allowed space region will be de- 
termined by the smaller of the two corresponding length 
scales — except for some narrow intermediate region 
when the two length scales are equal or similar. A simple 
as well as reasonable formula for estimating such a joint 
length scale from two ones is 



1 



I 2 

joint 
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I 2 



1 

I 2 

'2 



which in our case reads 
1 1 



1 



(2a B ) 2 (i?oA) 2 



(A5) 



(A6) 



Concerning the space region of the positron, we make the 
simplest choice to consider it unaltered with respect to 
the case of the free positronium, i.e., to be 2a B . The 
average distance between the electron and the positron, 
r op , also has to be estimated: For this we can use the 
formula 



(A7) 



as being a "half-way" between the minimal distance \r\ - 
r 2 1 and the maximal one r\ + r 2 : 



r\ - r 2 \ = \ r{ - 2r x r 2 +r^< 



< y r 2 + 2r x r 2 +r 2 =n+r 2 (A8) 

t 2 > 0). Putting all these together, we estimate the 
energy as 



1 



1 



^(approx) = Ke + R + V{ ) = 

2m G r G 2m c (2aB) 

„2 ^ 

(A9) 



47re ^(r c )2 + (2a B )2 

As can be seen on Fig. this formula provides a pretty 
good approximation for the ground-state energy of the 
confined positronium. 

We can apply the same style of approach for the case 
of the hydrogen atom in a hole as well. There, the joint 
length scale of the length scale a B of the free hydrogen 
and the one of the hole, Rq/tt, is given by 



(A10) 



In the potential energy this same length scale will ap- 
pear, since the proton stands in the center of the hole. 
Consequently, our estimate for the energy of a hydrogen 
atom in a hole is 



E, 



(approx) 

H 



n 2 l 



i 



(All) 

2m c rj| 47T£ r H 

Fig. [fj] shows that our reasoning gives a fairly good ap- 
proximation for this system as well. 
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APPENDIX B: THE PROBABILITY Pout FOR 
LARGE HOLES 

An explanation for the found large- Rq asymptotic be- 
havior p ou t ~ Rq 3 can be given by the following argu- 
ment: Let us rewrite Eq. © as 

oo Ro Xi 

Pout = 87r 2 J R 2 dR J r\ dr c J sin X d X |*(r e , R, x)| 2 , 

xi 

(Bl) 

where the integration boundaries Xi ^ and X2 < 
7r for x are determined from the condition r p = 
^/r^ + R 2 + 2r c Rcos\ > Ro- Expanding \^(r c , R, x)\ 2 
corresponding to Eq. (| 1 31) and also expanding all ips hi it 
in terms of <ps, let us examine the integral of one term, 

(r„) 4 X) (r e ) • V f (R) <pf (R) ■ <p™ ( X ) ^ (x) • (B2) 

Concerning the variable r e , the main contribution to the 
integral must come from r c i?o , within the range 
\r c — Rq\ ~ £ , since otherwise i? has to be around i?o to 
ensure r p > Ro , but that is exponentially suppressed by 
e -2-R/C w g-2fl /4 m ^ ne integrand. Now, for r c » i?o and 
i? <C i?oi the condition r p > i? is fulfilled by x in the 



region < x < 7r /2 • The integral of sinx^L (x) f i' ix) 
over this region gives an i? -independent numerical fac- 
tor. Turning to the r e -integration, for r c « R , the 
factor (^'(Ve) (r e ) behaves as 

rl(rl-Rly+ 1 ' = r e 2 (r e+J R ) i+i '(r e -Ji:o) i+i ' 

w i?2(2i? )^'(r c -i?o) J+l ' (B3) 

so its integral within the range < \r c — Rq\ ~ ^ is 
proportional to +2 ^ l+i +1 . The normalizing factors 
standing before (pj(r e ) and (p^}\r c ) bring in a factor 
p^~(2i+2i +3)^ regu ^ m g j n a behavior i? +1 \ The last 
integration, with respect to R, is an independent one, 

J R 2 dR ipf ] {R)ipf{R), (B4) 
o 

with a result depending only on £ so it does not modify 
the i?o-asymptotics (note that, for large Rq, £ is indepen- 
dent of the hole size, it tends to the constant 2ae). Hence, 
we can see that the strongest i?o-asymptotics among the 
terms in p out is caused by the i = %' = 1 term, and is 
found to be R^. 3 ■ 
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